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Abstract 

This article is concerned with the numerical simulations of perfect crystals. We study 
the rate of convergence of the reduced Hartree-Fock (rHF) model in a supercell towards the 
periodic rHF model in the whole space. We prove that, whenever the crystal is an insulator 
or a semi-conductor, the supercell energy per unit cell converges exponentially fast towards 
the periodic rHF energy per unit cell, with respect to the size of the supercell. 


1 Introduction 

The numerical simulation of the electronic structure of crystals is a very active research area in 
solid state physics, materials science and nano-electronics. When the crystal is perfect, a good 
approximation of its electronic ground state density can be obtained by solving a mean-field non¬ 
linear periodic model set on the whole space. Using the Bloch transform [14, Chapter XIII], we 
can recast such a problem as a continuous family of compact problems indexed by points of the 
Brillouin-zone. In practice, the compact problems are solved on a discretization of the Brillouin- 
zone. There is therefore an inherent error coming from the fact that the Brillouin-zone is sampled, 
and it is not obvious a priori whether this error is small, due to the nonlinearity of the problem. 
It has been observed numerically since the work of Monkhorst and Pack [12] that this error is 
indeed very small when the discretization is uniform, and when the crystal is an insulator or a 
semiconductor. To our knowledge, no rigorous proof of this fact was ever given. This article 
aims at proving why it is indeed the case in the reduced Hartree-Fock (rHF) model, which is a 
Hartree-Fock model in which the exchange term is neglected. This model was studied in [4, 5]. 

A crystal is modeled by a periodic nuclear charge distribution p-pev The corresponding rHF 
energy per unit cell is denoted by Ipi”. When numerical calculations are performed over a regular 
discretization of the Brillouin-zone, this amounts to calculate the energy on a supercell, i.e. on a 
large box containing L times the periodicity of piper hr each direction (for a total of unit cells 
in the supercell), and with periodic boundary conditions. The rHF energy on a supercell of size L 
is denoted by so that the corresponding energy per unit cell is L 

It was proved in [4] that converges to Ipe” as L goes to infinity, when the crystal is 

an insulator or a semiconductor. However, following the proof in [4] , we find a rate of convergence 
of order L~^, which is well below what is numerically observed. Our main result is that, if the 
crystal is an insulator or a semiconductor, then there exist constants C € K+ and a > 0, such 
that 

VL e N*, < Ce-°‘^. (1.1) 

We also prove that the supercell electronic density converges exponentially fast to the periodic 
rHF electronic density, in the norm. To prove such rates of convergence, we recast the 

problem into the difference between an integral and a corresponding Riemann sum, and show 
that the integrand is both periodic and analytic on a complex strip. Similar tools were used 
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in [7, 8, 10, 2, 13] to prove that the Wannier functions of insulators are exponentially localized. 


This article is organized as follows. In Section 2, we recall how the rHF model is derived, 
and present the main results. In Section 4, we apply the Bloch theory for both periodic models 
and supercell models. The proofs of the main results are postponed until Section 5. Finally, we 
illustrate our theoretical results with numerical simulations in Section 7. 

Throughout this article, we will give explicit values of the constants appearing in the inequal¬ 
ities. These values are very crude, but allows one to see how these constants depend on the 
parameters of the electronic problem. 

2 Presentation of the models 

A perfect crystal is a periodic arrangement of atoms. Both the nuclear charge density /iper and 
the electronic density are 7^-periodic functions, where 77 is a discrete periodic lattice of Let F 
be the Wigner-Seitz cell of the lattice, and let F* be the Wigner-Seitz cell of the dual lattice 77* 
(one can also take F* as the first Brillouin-zone). For instance, for 77 = aZ^, F = [—a/2,a/2)^, 
77* = {2tt! a)1? and F* = [—tt/o, tt/o)^. For R G 77, we let tr be the translation operator 
on defined by (tr/)(x) := /(x— R). 

We will assume throughout the paper that the nuclear charge density /Xper is in Lpgj.(F) for 
simplicity, but distributions with singularity points may also be handled [1]. 

2.1 The supercell rHF model 

In a supercell model, the system is confined to a box F/, := LT with periodic boundary conditions. 
We denote by l) the Hilbert space of locally square integrable functions that are L77- 

periodic. The Fourier coefficients of a function / e Lper(r l) are defined by 

Vk e L-^77*, c^(/) = f /(x)e"'*" ’‘dx, 

d l\ Jfl 

so that, for any / S Lp^j.{T l), 

/W = a.e. and in L2^,(ri,). 


The set of admissible electronic states for the supercell model is 

Vl ■■= e S{LI^^(Tl)), 0 < 7l < I, TrL 2 ^ 7 r^) (yz,) Trij^yr^,) (-AlJl) < ooj , 

where S{'H) denotes the space of the bounded self-adjoint operators on the Hilbert space H. Here, 
Trz^g^riri,) {—AlJl) is a shorthand notation for 


3 

TrL2_^7ri) (-Al7l) := ^ Tri^^yr^) {Pj,LlLPj,L ), (2.1) 

where, for 1 < j < 3, Pj^L is the self-adjoint operator on Lpgj.(r l) defined by cl^{Pj,Lf) = kjC^{f) 
for all k = {ki,k 2 ,k^) G L~^TZ*. Note that Ck(—Az/) = |kpCk(/) for all k G L~^TZ*. 

We introduce the L77-periodic Green kernel Gl of the Poisson interaction [11], solution of 

/ —AGl = 47r (X)ke7?, — l) 

1 Gz is L77-periodic. 
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The expression of Gl is given in the Fourier basis by 


A-tt _ pikx 

Gnx) = c, + — y: (2.2) 

' kGL-i7?,*\{0} ' ' 

where cl = fp^ Gl- The constant cl can be any fixed constant a priori. In one of the 

first article on the topic [11], the authors chose to set cl = 0, but other choices are equally valid 
(see [4] for instance). This is due to the fact that cl does not play any role for neutral systems. 
We choose to set cl = 0 for simplicity. The supercell Coulomb energy is defined by 

V/,5 e Tper(rL), T>L(/,g) :=^(/*ri. GL)(x)5r(x)dx, (2.3) 

where (/ Gl){x) := /(y)Gi(x — y)dy. We recall that the map p p *ri. Gp is continuous 
from L^er(rL) to L^eri^L)- 

Any € Vl is locally trace-class, and can be associated an L72.-periodic density p^^ S 
Lpgr(ri). For e Vli the supercell reduced Hartree-Fock energy is 

■S’l'’"(7l) := ^Tri2^^(r^)(-Ai7i) -b - /iper,P7i, “ Mper)- (2.4) 

The first term of (2.4) corresponds to the supercell kinetic energy, and the second term represents 
the supercell Coulomb energy. The ground state energy of the system is given by the minimization 
problem 

= inf |£'£-'(7l), 1l G Vl. ^ ^ Mperj • (2.5) 

Using techniques similar to [4, Theorem 4], the following result holds (we do not prove it, for the 
arguments are similar to the ones in [4]). 

Theorem 2.1 (Existence of a supercell minimizer). For all L S N*, the minimization prob¬ 
lem (2.5) admits minimizers. One of these minimizers 7l,o satisfies tu'jl.o = Il.oTr. All mini- 
mizers share the same density p-y^ g, which is TZ-periodic. Finally, ^l.o satisfies the self-consistent 
equation 

{ 7l,o = 1 {Hl,o < £f) + ^ 

I Hlp =-\^l + Vlp ( 2 . 6 ) 

[_ = {P'Tl.o ~ Mper) *r Gl. 

where Hl.o acts on Lpgj.(ri) and 0 < <5 < 1{F[lp = Sp) is a finite rank operator. 

Here, Sp is the Fermi level of the supercell model. It is chosen so that the charge constraint 
in (2.5) is satisfied. 

Remark 2.2. The LTZ-periodic density of the minimizers g is actually TZ-periodic. It is unclear 
that such a property should hold for more complex models (e.g. Kohn-Sham models). This is the 
reason why we state our results for the rHF model. We believe however that similar results should 
hold true for more complex systems, provided that the supercell density is TZ-periodic for each size 
of the supercell. 

2.2 The reduced Hartree-Fock model for perfect crystals 

The rHF model for perfect crystals, or periodic rHF, has been rigorously derived from the rHF 
model for finite molecular systems by means of a thermodynamic limit procedure by Catto, Le 
Bris and Lions [5]. In [4], Cances, Deleurence and Lewin proved that the same periodic rHF model 
is the limit of the rHF supercell model as the size of the supercell goes to infinity. 


3 



We introduce the set of admissible density matrices 

Vper ■■= {l G 0<7<1, VRe7^, tr7 = 7Tr, Tr ( 7 ) + W (-A 7 ) < 00 } , (2.7) 

where Ty denotes the trace per unit volume. For any locally trace class operator A that commutes 
with 77.-translations, it reads 

W (71) := lim 4Tr(lLrAlLr). (2.8) 

L-s-oo 

The trace per unit volume Tp can also be defined via the Bloch transform (see Equation (4.6) 
below). Here, Tp(—A 7 ) is a shorthand notation for 

3 

W(-A 7 ) :=^W(H, 7 P,), 
i=i 

where Pj = is the momentum operator in the direction. The Coulomb energy per unit 

volume is defined by 

V/,ff e Lper(r), Di{f,g) := J^{f *rGi){x)g{yi)dx, (2.9) 

where Gi was introduced in ( 2 . 2 ). 

Any 7 G Vper is locally trace-class, and can be associated an 77.-periodic density p-y G Lpg,.(r). 
For 7 G Pper, the reduced Hartree-Fock energy is given by 

^pep(l) ~ ~ A^per) ' ( 2 - 10 ) 

The first term of (2.10) corresponds to the kinetic energy per unit volume, and the second term 
represents the Coulomb energy per unit volume. Finally, the periodic rHF ground state energy is 
given by the minimization problem 

IpiP ■= inf |f^eT(7), 7 G 'Pper, j^P'r= Mperj ■ (2.11) 

It has been proved in [4] that the minimization problem (2.11) admits a unique minimizer 70 , 
which is the solution of the self-consistent equation 

(70 = l{Ho < ep) +S 

I Ho = -iA + Eo (2.12) 

^0 (P70 Pper) *r Gi^ 

where Hq acts on L^(IR^) and 0 < (5 < 1{Hq = ep) is a finite rank operator. Here, the Fermi 
energy sp is the Lagrange multiplier corresponding to the charge constraint /p pyg = /p Pper- We 
make the following assumption: 


(Al) The system is an insulator, in the sense that Hq has a spectral gap around sp.\ 

In particular, 5 = 0. 

3 Main results 

Our main results are concerned with the rate of convergence of supercell models towards corre¬ 
sponding periodic models. We first prove the exponential rate of convergence in a linear setting, 
where the mean-filed potential E is a fixed 7^-periodic function: E G L“j.(r). We then extend our 
result to the nonlinear rHF model, where the external potential is the solution of the self-consistent 
equation (2.6) or (2.12). 

We start with the linear case. The proof of the following proposition is given in Section 5. 
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Proposition 3.1 (Convergence rate of the linear supercell model). Let V G L“j.(r) be such that 
the operator H = —^A + V acting on has a gap of size g > 0 centered around the Fermi 

level Ef- Then, for any L G N*, the operator Fl^ := —+ V acting on Lpgj.(ri) has a gap of 
size at least g around ep ■ Let 

'y = t{H < Ef) and 7 l = 1 < Ef) ■ (3.1) 

Then, 7 G Vper and 7 ^ G Vl, and there exist constants C G K"*" and a > 0, that depend on the 
lattice TZ, ||y||icxj, g and Ef only, such that 

VL G N*, |Tl ( 7 ^?) — Tr j- | < Ce““^ (ground state energy per unit volume) (3.2) 

and 

VL G N*, \\p^ — < Ce~°‘^ (ground state density). (3.3) 

In a second step, we will use the projectors 7 and 7 ^ obtained for well chosen potentials V 
as candidates for the minimization problems (2.11) and (2.5) respectively. We have the following 
result (see Section 5.5 for the proof). 

Corollary 3.2. With the same notation as in Proposition 3.1, there exist constants C G K’*' and 
a > 0, that depend on the lattice TZ, |jC||i.x,, g and Ep only, such that 

VL e |£:p^|r(7) - L-3£:^-="(7i)| < Ce-“^. (3.4) 

We are now able to state our main result for the rHF model. The proof of the following theorem 
is given in Section 6 . In the sequel, we denote by B{E) the set of bounded operators acting on the 
Banach space E. 

Theorem 3.3 (Convergence rate of the rHF supercell model). Under hypothesis (Al), there 
exist C G K’*' and a > 0 independent of L such that the following estimates hold true: 

• convergence of the ground state energy per unit volume: 

VL e N*, - 7^1” I < Ce-“^; 

• convergence of the ground state density: 

VLeN*, ||p^,,„-p^„||L~(r) <Ce-“^; 

• convergence of the mean-field Hamiltonian: 

VL e N*, ||i?L — H‘oI!b(l2(b3)) < 

where Hp := — iA+(p.yj^ ^ — Ppei) *rGi and H := —}^A + {p-^^ — /rper)*rGi are acting on L^(IR^). 

The fact that the supercell quantities converge to the corresponding quantities of the periodic 
rHF model was already proved in [4, Theorem 4]. However, following the proof of the latter article, 
we only find a O (L“^) convergence rate. 

The proof of Proposition 3.1 and Theorem 3.3 rely on Bloch transforms. 


4 Bloch transform and supercell Bloch transform 

4.1 Bloch transform from to Lpgj,(r*, L^(r)) 

We recall in this section the basic properties of the usual Bloch transform [14, Chapter XHI]). Let 
(ai, a 2 , aa) be a basis of the lattice TZ, so that TZ = Zai + Za 2 + Zaa. We define the dual lattice 
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TZ* by TZ* = Za* + Za^ + Zag where the vectors a* are such that a* • a^ = 2TT5ij. The unit cell 
and the reciprocal unit cell are respectively defined by 

r := {oiai + a2a2 + Oaag, (oi, a2j c^a) € [~l/2,1/2)^} , 

and 

r* := {aia* + aaaj + aga*, (oi, 02 , aa) G [-1/2,1/2)3} ^ 

Note that T* differs from the first Brillouin-zone when the crystal is not cubic. We consider the 
Hilbert space L^(r*, Lpgj.(r)), endowed with the normalized inner product 


(/(q,x),g(q,x))i2(r.,i2^^(r)) := ^ J /(q, x)g(q, x) dxdq. 


The Bloch transform is defined by 

Z: L^R^) ^ (D) 


w I— 


(Zu;)(q,x) :=u;q(x) ^ e-‘'i'(^+*^)zu(x + R). 


(4.1) 




Its inverse is given by 

z-W L^{r*,L%,{r)) ^ 

u;q(x) !->• (Z“^w)(x) :=-f e‘'*'^ri;q(x) dq. 


It holds that Z is an isometry, namely 


ll^^lli2(ryL2„(r)) ^ l(2^w^)(q, x)|^ dx dq = || ^lli2(R3)- 

For m gTZ* , we introduce the unitary operator Um acting on Lpgj.(r) defined by 

Vme7^^ yf€Ll,,(r), (C7^/)(x) = e-™V(x). (4.2) 

From (4.1), it is natural to consider Zw as a function of (K^^ Lpgj.(r)) such that 

yw€L‘^{M.^), Vme7^*, VqSF*, (Zw) (q + m, •) = Wq+m = tlmiCq = Um (Zu;(q, •)). 

(4.3) 

Let A with domain T>{A) be a possibly unbounded operator acting on Lpg,.(r). We say that A 
commutes with 7^-translations if trH = Htr for all R G 7^. If H commutes with 7?,-translations, 
then it admits a Bloch decomposition. The operator ZHZ“^ is block diagonal, which means 
that there exists a family of operators ( 4 lq)qgp. acting on Lpg,.(r), such that, if / G L‘^(U.^) and 
g G 'D(A) are such that / = Ag^ then, for almost any q G F*, g^ G Lpgj.(F) is in the domain of 
Aq, and 

/q = ^q5q- (4-4) 

In this case, we write 


ZHZ“^ = -r Hqdq (Bloch decomposition of A). 


From (4.3), we extend the definition of Hq, initially defined for q G F*, to q G with 

VmG7^^ VqGF*, A^+^ = U^A^U-\ (4.5) 

so that (4.4) holds for almost any q G If H is locally trace-class, then A^ is trace-class on 
Lper(r) for almost any q G xhe operator A can be associated a density pA, which is an 
7^-periodic function, given by 
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where pA^ is the density of the trace-class operator A^. The trace per unit volume of A (defined 
in ( 2 . 8 )) is also equal to 


It(^) = Tri2^^(r) (v4q) dq. 


(4.6) 


4.2 Bloch transform from Lp^^{Tl) to Lp^^{T)) 

We present in this section the “supercell” Bloch transform. This transformation goes from Lpgj.(ri) 
to £^(Al, Lpgj.(r)), where := n T*, i.e. 


At, := 


2ki * 2fc2 * ^ 2^3 ^ 

— —a.2-\ — \kijk2,k^) & 


— L + T] —L + 7] 


-1 


(4.7) 

with 77 = 1 if L is odd, and 77 = 0 if L is even, so that there are exactly points in A^. Similarly, 
we define TZl ■= 77.(7 T l, which contains points of the lattice 77. The supercell Bloch transform 
has properties similar to those of the standard Bloch transform, the main difference being that 
there are only a finite number of fibers. We introduce the Hilbert space ^^(Al, Lpgj.(r)) endowed 
with the normalized inner product 


(/(Q, x), g{Q, x))f2(A^_i2^_^(r)) 


1 

lA 


/ /(Q>x)g(Q,x)dx. 

QeAi, 


The supercell Bloch transform is defined by 

Za: ^ eHAL,Ll,Ar)) 

w H> (Zlw)(Q, x) := 77;q(x) := ^ e“'‘^'^’‘’*'^^7(;(x -1-R). 


Its inverse is given by 

Z-i; e{AL,LlA^)) ^ 

wq{x) ^ ^ 

QeAi 

It holds that is an isometry, i.e. 

QeAi 

We can extend Z to £°° (L“^77*, Lpgr(r)) with 

Vtu e Lpg,(rL), Vm e 77*, VQ G Al, wq+m = UmWQ, 
where the operator t/m was defined in (4.2). 

Let A^ with domain V (A^) be an operator acting on Lpgj.(rA). If A commutes with 77- 
translations, then it admits a supercell Bloch decomposition. The operator Z^A^Z^ is block 
diagonal, which means that there exists a family of operators (AQ)QgAi acting on Lpg^(r) such 
that if / = A^g with / G Lpgj.(rA) and g G 'D{A^), then for all Q G A^, 


/q = ^q5q- 


(4.8) 


We write 


ZlA^zj;^ 


Aq (supercell Bloch decomposition of A^). 

QeAi 
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(4.9) 


The spectrum of can be deduced from the spectra of (^Q)Qg^ with 

= U (^q) • 

QeAi, 

Similarly to (4.5), we extend the definition of Aq to L~^TZ* with 

VmG7^^ Vqer, AQ+^ = UmAQU-\ 

so that (4.8) holds for all Q G L~^TZ*. 

Finally, if the operator A^ is trace-class, we define the trace per unit volume by 


(4.10) 


QeAi 


and the associated density is given by ^ X] where is the density of the 


trace-class operator Aq. 


5 Proof of Proposition 3.1: the linear case 

The proofs of Proposition 3.1 and Theorem 3.3 are based on reformulating the problem using the 
Bloch transforms. Comparing quantities belonging to the whole space model on the one hand, and 
to the supercell model on the other hand amounts to comparing integrals with Riemann sums. 
The exponential convergence then relies on two arguments: quantities of interest are 7^*-periodic 
and have analytic continuations on a complex strip, and the Riemann sums for such functions 
converge exponentially fast to the corresponding integrals. 

We prove in this section the exponential convergence of Proposition 3.1. 


5.1 Convergence of Riemann sums 

We recall the following classical lemma. For A > 0, we denote by 

5^ := {zGC^, |Im(z)U < A} + 

If i? is a Banach space and d G N*, an Fl-valued function F : If C — >■ i? is said to be (strongly) 
analytic if (VzF)(z) exists in E'^ for all z G fl. In the sequel, we assume without loss of generality 
that the vectors spanning the lattice TZ* are ordered in such a way that |a* I < |a^| < |a5|. 

Lemma 5.1. Let / : —>■ C be an TZ* -periodic function that admits an analytic continuation on 

Sa for some A > 0. Then, there exists C G K’*' and a > 0 such that 


VL G N*, 


/(q)dq- 7^ X 


QeAi, 


< Co sup |/(z)|e 
zGSa 


The constants may be chosen equal to 


a = (2/3)7rA|a3| and Cq = 2 


3 -b e 


-2a 


(1-e-^r 


(5.1) 


Proof of Lemma 5.1. Let cr(/) := /p. /(q)e '^ "^dq be the Fourier coefficients of /, so that 

/(q) = X 

rgt?. 
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It holds 


r* 


/(q)dq-^ /(Q) 

QeAi, 


Y E CR(/)e‘Q 


R 


QgAl Re7^ 


E I E 


Q R 


RG7?,\{0} 


QeAi 


By noticing that 


we obtain 


QgAx, 


0 if R^LTZ 
otherwise , 


/ /(q)dq-y^ Y /(Q) 

= 

E Clr(/) 

QeAi 


Re7?,\{0} 


(5.2) 


If / is analytic on Sa, we deduce from /(q) = X^rgk analytic continuation 
of / is given by 


-Ry 


Vq e KZ Vy e [-A, /(q + iy) = ^ CR(/)d^"i e 

so that {cR(/)e“*^'^}j^g.j^ are the Fourier coefficients of the 7?.*-periodic function q n- /(q + iy)- 
In particular, 

VReT^, VyG[-A,A]^, |cr(/)| < sup |/(q + iy)|e*^'5". (5.3) 

qer* 

We make the following choice for y. We write R = fciai + k 2 a 2 + ^ 3^3 with ki, k 2 , k^ G Z, and we 
let I < TO < 3 be the index such that \km\ = \kj\^- Choosing y = —sgn(kin)A|a 3 |“^a)(j G [—A,A]^, 
leads to 

|cr(/)| < sup |/(z)|e-2-Aa3|-bfeL < sup |/(z)| e-“l'=di 
zGSa zGSa 

where we used the inequality |k|oo > (l/3)|k|i, and we set a = (2/3)7rA|a3|“^. Note that the 
Fourier coefficients of / are exponentially decreasing. We conclude with (5.2) and the inequality 

g-QL|kq 


E CiR(/) 

RG7?,\{0} 


< E E 

RG7?,\{0} kGZ3\{0} 


< sup |/(z) 
zGSa 


^g-aL|fe| 




= »p |/(.)| e- 

2 &Sa \ (I — e “)^ / 


□ 


5.2 Analyticity and basic estimates 

The exponential rates of convergence observed in (3.2), (3.3) and (3.4) will come from Lemma 5.1 
for appropriate choices of functions /. In order to construct such functions, we notice that H and 
defined in Proposition 3.1 commute with 7?.-translations, thus admit Bloch decompositions. 
From 


VqGMZ (-A)q = |-iVi + q|' = 


3 

E 

i=i 


{Pj,i + < 13 ) and VQ G L-^n* 


(-Al)q = (-Ai) 


Q 
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where Vi denotes the gradient on the space Lpgj.(r) and Ai was defined in ( 2 . 1 ), we obtain 

ZHZ~^ = J Hqdq with i/q |-iVi + + = - (-Ai - 2iq • Vi + |qp)+y, (5.4) 

and 

QeAi, 

In other words, for all Q in Al, = Hq. In addition, the spectrum of H can be recovered 

from the spectra of (^^q)qgr* [14, Chapter XIII] 

<^iH) = IJ cr(7?q) 

qer* 

Together with (4.9) we deduce that, since H has a gap of size g centered around ep, then has 
a gap of size at least g around sp¬ 
in the sequel, we introduce, for z G C^, the operator (we denote by := lor simplicity) 

iJz := i (-Al - 2iz • Vi + z^) + 14 acting on Lpg,.(r). (5.5) 

With the terminology of [9, Chapter VII], the map z is an holomorphic family of type (4l). 

Let S := inf a{H) be the bottom of the spectrum of H. We consider the positively oriented simple 
closed loop U *^2 U *^3 U “^4 in the complex plane, consisting of the following line segments: 

= [ef — i, ef + i], ^2 = [ef +i, X — 1 + i], = [S — 1 + i, S — 1 — i] and ^4 = [E — 1 — i, — i]- 


'^2 


^3 


7 

^1 a{H) 


5 

L 

ef 


Figure 1: The loop ^. 


The projectors defined in (3.1) can be written, using the Cauchy’s residue theorem, as 

1 jf dA , 1 r dA 

Together with (5.4), it follows that 7 and jp commutes with 7?.-translations, with 
2 'yZ~^ = J 7 qdq and ZL"fL2l^ = — 0 7Q, 


QeAi 


where we set 


Vq G 


7q := 


dA 


2 iF X- Ha 


For Q G TZ*, it holds (7l)q = 7q- The analytic continuation of (5.7) is formally 


VzG^, 1.-.= — 


dA 


1 


2iF jAg X- H„ 


(5.6) 


(5.7) 
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The fact that A — H^, is indeed invertible, at least for z in some Sa for A > 0 is proved in the 
following lemma. For z S C^, and X we introduce 


gi(A,z) := (1 - and g 2 (A,z) := ^ (1 - Ai). (5.8) 

Lemma 5.2. For all q G and all A € the operator A — iJq is invertible, and there exists a 
constant Ci G K“'" such that, 

Vqer, VAe-r, ||i?i(A,q)||g(^.^^(r)) <Ci and ||i? 2 (A, q)||e(i.^^(r)) < Ci. (5.9) 

Denoting by |r *|2 := sup{|q| 2 , q € F*}, we can choose 

c.=4 + ^ + 4|rii + 8||r|U-+8.p 
mm(l, 5 ) 

Moreover, there exists ^4 > 0 such that, for all z G Sa and all X G ^, the operator X — Fl^ is 
invertible, and there exists a constant C 2 G such that 


Vz G F* + i [—A, A]^, VAg^, ||i3i(A, < (72 and ||B 2 (A, z)||g^p 2 ^^^pjj < ( 72 . 

(5.11) 

IFe can choose 

A = 2Ci(1 + |F*|2) ) ^2 = 2Ci. (5.12) 

This lemma was proved in the one-dimensional case by Kohn in [10], and similar results were 
discussed by Des Cloizeaux in [7, 8 ] . 

Remark 5.3. The bounds (5.9) and (5.11) are not uniform for q G (only for q G F*/ This 
comes from the fact that, for m G TZ*, 


1-Ai 

^ / e ^ 

1-Ai 

l-f|m|^ 

1 - 1 - (-iVi - 1 - m)2 


1 - 1 - (-iVi - 1 - m)2 

|F|i /2 / |r| 


Remark 5.4. From Lemma 5.2, we deduce that {lz)y,^SA analytic family of bounded oper¬ 

ators. Since 7 q is an orthogonal projector for q G i.e. jq = 7 q 7 q, we deduce that 7 z = 7 z 7 z 
for all z G Sa, so that 7 z is a (not necessarily orthogonal) projector. Also, Tr( 7 z) is a constant 
independent of zG Sa- 

Proof of Lemma 5.2. From the inequality jap < 2|a -I- -I- 2|&p, we get that, for q G F*, it holds 


|—iVi -I- q|^ -I- jqp > — 5 A 1 . We deduce that 

VqGF*, Hq>-\/X,-\\T*\l-\\V\\L^. (5.13) 

We first consider the part of the contour ^ (see Figure 1). It holds 

VAg^i, VqGF*, \Hq-X\^>\Re{Hq-Xt = \Hq-eF\\ (5.14) 

Since |i7q — e_F| > g/2, we get 

VAg-^i, VqGF*, |i7q-A|>g/2. (5.15) 

On the other hand, from (5.13) and (5.14), it holds that 

VAg^\ VqGF*, \Hq-X\>Hq-eF>-\Ax-^\T*\l-\\V\\L^-eF. (5.16) 
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Combining (5.15) and (5.16) leads to 


VM>0, VAe-^i, Vqer*, (M + 4)|i/q-A| >-Ai+M|-2|r*|2-4||C||Loo 
Choosing M = (2 + 4|r*|| + 8||1A||lc=o + 8£F)/g gives 


VAG‘r\ qer*, 




A ^ 2 + 4|r*|g,+8||y|U^+8£F 

“V 9 


(1 - Al), 


which proves (5.9) for A G '^i. The inequalities on the other parts of ^ are proved similarly, the 
inequalities (5.15) and (5.16) being respectively replaced by their equivalent 

VA G ^2 U <^ 4 , l-ffq - A|2 > |Im (i?q - A)|^ > 1 and |i?q - A| > i/q - S - 1 > i/q - ep, 

VAg^ 3 , |i?q-Ap > |Re (i/q-A)|^ > 1 and |i?q - A| > i/q - E - 1 > i/q - ep. 

This proves (5.9). We now prove (5.11). For z = q + iy G with q G T* and y G one can 
rewrite (5.5) as 


= i?q + y • Vi + iq - y - -|yp = i^q + y • Vi - -y + iq . 


In particular, 


A - Ff, = A - Ffq + Ffq - FT, = (A - Ffq) 1 - (A - i/q) 


y 1 Vi - -y + iq 


= (A-i7q) (^l-il2(A,q)^—^ 

For |y|oo < 1, we have 

< |ylc 


y • (- 2 ^+ 


(5.17) 


1-A, 


y ■ ( Vi - -y + iq 


|Vi| 


1-Ai 


+ 2 !^!°° ■*" 1^1°° 1 - l^loo (1 + |F*|2) ■ 


Together with (5.9), we obtain that for all |y|oo < A := min (l, (2Ci(l + |r*| 2 )) ^), 


-82 (A, q) 


1-Ai 


y ■ 1 Vi - -y + iq 


1 

< . 

- 2 


As a result, from (5.17), we get that for all q G T* and all y G [—A, A], the operator A — H^. is 
invertible, with 


X-H, 


(1-Ai) 


< C 2 ■= 2Ci and 


ecHi) 


(1-Ai) 


A-i7. 


< C 2 . 


B(Wi) 


□ 

For z G Sa, introduce the operators Bi{z) and B 2 {z) respectively defined by 

Bi(z) := (1 - Ai) 7 z and ^ 2 ( 2 ) := 7 z(l - Ai). (5.18) 

In the sequel, for k G N*, we denote by 6 fc('H) the fc-th Schatten class [15] of the Hilbert 
space H ; 61 (H) is the set of trace-class operators, and 62 (H) is the set of Hilbert-Schmidt 
operators. From Lemma 5.2, we obtain the following result. 
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Lemma 5.5. There exists a constant C 3 € such that 


VzgT* +i[-A,A]^, Bi{z) 
The value of C 3 can he chosen equal to 


r(r)) 


< C 3 and 


B 2 {z) 




<^3. 


C3 = -Ci{3 + eF + \\V\\L^). 


Also, for all z G F* + i[—A, yl]^, the operator 7 ^ is trace-class, and 


llT'zllsi(L 2 _^,(r)) — ^4 with C 4 — Cl ^ fi + lkp') ■ (5.19) 

kGT?,* ^ 11/ 

Proof. The first assertion comes from the fact that 

.Bi(z) = Bi{\,z) dA, 

and the fact that \‘t^\ = 6 2{£f — E) (see Figure 1). Note that since |—iVi + q|^ > 0, it holds 

S > — To get the second assertion, we note that 7 z is a projector, so that 

7z = 7z7z = B2{z) 

The operator (1 — Ai)“^ being trace-class, with 


k|2 


1(1-A)- 


lei(L2 (r)) 


= E 

kGT?,* 


1 + 



we obtain (5.19). □ 

5.3 Convergence of the kinetic energy per unit volume 

The kinetic energy per unit volume of the states 7 and defined in (3.1) are respectively given 
by 

iFper := dr (-A 7 ) and Kl ■= Tcl ■ 

Using the Bloch decomposition of 7 and 71 , in (5.6)-(5.7), and the properties (4.6) and (4.10), we 
obtain that 

3 3 

iFP- = ^ / iF,(q) dq and = ^^(Q) 

i=l j = l qSAz, 

where, for 1 < j < 3, we introduced the function 

iFj : q 9 H> Tr^^^^p) {{Pj + {Pj + qj)) ■ 

Here, we denoted by Pj := Pij for simplicity. Recall that the operator Pij was defined in (2.1). 
The error on the kinetic energy per unit volume FfP®"' — is therefore equal to the difference 
between integrals and corresponding Riemann sums. In the sequel, we introduce, for 1 < j < 3, 
the function 

Vz e Sa, K,{z) := Tri 2 ^ 7 r) {fPj P zj) 7z (Pj + zj)). 
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Lemma 5.6 (Exponential convergence of the kinetic energy). For all 1 < j < 3, the function Kj 
is TZ-periodic, and admits an analytic continuation on Sa, where A > 0 was defined in (5.12). 
Moreover, it holds 


sup |itrj(z)| < Cs where C 5 = Mr *|2 + ^ + ^ ) ClC^. (5.20) 

As a consequence, from Lemma 5.1, it holds 

< CoC5e-°‘^, 

where Cq € K’*' and a > 0 were defined in (5.1). 

Proof. The 7^-periodicity comes from the covariant identity (4.5). To prove the analyticity, it is 
enough to prove that {{Pj + Zj)'^z.[Pj + Zj)) is a trace-class operator for all z S T* -|-i[—A, 

We only consider the case j = 1 and fc = 1, the other cases being similar. We have 

dzi {{Pi + ^i)7z(^i + ^ 1 )) = 7z(^i + ^ 1 ) + (^1 + ^i)(52i7z)(^i + zi) + {Pi + Zi)'jz- (5-21) 

We first show that (Pi -I- ^ 1)72 is a bounded operator. We have 

(Pi -I- 2:i)7z = 

where Pi was defined in (5.18). From Lemma (5.5) and the fact that (Pi -I- zi){l — Ai)“^ is a 
bounded operator, we deduce that (Pi -I- 2 : 1)72 is bounded. The proof is similar for the operator 
7 z(Pi + Zi). We now turn to the middle term of (5.21). Since 7 ^ is a projector, it holds 7 ^ = 7 z 7 z- 
We obtain 


(Pi -k 2 :i)(i9zi7z)(Pi + zi) = {Pi + zi) [7z(9^i7z) -f (9zi7z)7z] (Pi + zi) 

= [(Pi -I- zi)jz] 7z(9zi7z)(Pi + zi) + {Pi + 2:i)(5z,7z)7z [7z(Pi + ^i)] ■ 


We already proved that the operators (Pi -I- 21)72 and 7z(Pi + 21) were bounded. Also, 7z is a 
trace-class operator. To prove that (Pi -I-21) (9zi) 7z(Pi +21) is trace class, it is therefore sufficient 
to show that (Pi -I- Zi){dz,,jz) is bounded. We have 


(Pi -k2i)(a„,7z) 




X-H, 


dX 


1 

2i7r 


/ Pi + Zi 


Pi (A, 



2 

dA, 


which is a bounded operator. We conclude that dz^ {{Pi + 21 ) 72 (Pi + 21 )) is a trace-class opera¬ 
tor. Finally, for 1 < j < 3, Kj is an analytic function on Sa- 


To get the bound (5.20), we write that 


^j{'^) ~ {{^j ffi ■^i)7z(Pj + Zj)) — Tri;,2 (p) {{Pj + 2j)7z7z7z(Pj + Zj)) 


= TrL2 (D 


^^+^^Pi(z)72P2(z)-^^ 


" 1 -Ai 

The bound (5.20) easily follows from Lemma 5.5 and the estimate 


Vz G r -ki[-A,A]3, 


Pj P Zj 


1-Ai 


B(iSer(r)) 


^ Fj \ + 


1-Ai 


B(iSer(r)) 


< |r *|2 + Ai + v. 


□ 
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5.4 Convergence of the ground state density 

We now prove (3.3). The densities of 7 and 7 ^ defined in (5.6)-(5.7) are respectively 


P 7 := r P 7 q dq and '■= Jz P'fQ' 

QSAi, 

In particular, if IT is a regular 77,-periodic trial function, it holds that 

■= i = i ( 7 qW^) dq and := f ^ ^ Tr^^^^ ( 7 qIT) , 

so that the error is again the difference between an integral and a corresponding 

Riemann sum. We introduce, for W G Lpg^(r) the function 

Vz e Sa, Mw(2) := Tri 2 ^^ ( 7 ,W). (5.22) 

Lemma 5.7. For all W G Lpg,.(r), Mw defined in (5.22) is well-defined TZ*-periodic analytic 
function on Sa, where 4 > 0 was defined in (5.12), and it holds that 

2 


sup \Mwiz)\ < CelllTllLi^^jr) with Cq = Cj ^ 
zGSa 


ken- 


1 


l + |k|^ 


As a consequence, from Lemma 5.1, it holds that 

IIP7-P7i.llL~(r) ^CoCee-"-^, 
where Cq and a were defined in (5.1). 


(5.23) 


(5.24) 


Proof of lemma 5.1. We first prove that My/ is well defined whenever IT £ Lpg,.(r). For IT G 
Lpgr(r), we have 

Mw{z) = Tri,2^7r) {izW) = Tr^^^^r) (7zH"7z) = Tri2^dr) (i^(z)(l - Ai)-ilT(l - Ai)-ii?;(z)) . 
According to the Kato-Seiler-Simon inequality [15, Theorem 4.1]^, the operator (1 — Ai)“^\/|IT| 

2 \ 1/2 


is Hilbert-Schmidt {i.e. in the Schatten space 62 (Lper(r))), and satisfies 


(l-Ai)-!^^ 


e2(Aj„(r)) 




\keR* 


1 


l + |k|2 


11/2 


liw^llinr)- 


llh^ll 


Li(r), 


(5.25) 


It follows that (1 — Ai) ilT(l — Ai) 1 is in 6 i(Lpgj.(r)) with 

||(l-Ai) ilT(l-Ai) lei(L 2 „(r)) - (1 + [tp) 

The proof of (5.23) then follows from Lemma 5.5. 

Let us now prove that, for IT G ^^(r), My/ is analytic on Sa- To do so, it is sufficient to show 
that, for 1 < A: < 3, ^^^(yzlTyz) is a trace class operator. We do the proof for A: = 1. We have 

'92i(7zbF7z) = (92 i7z)IT7z +7zIT(5^i7z) 


+ 


1 


(1) (J) B 2 (A,z) 


^ -IT^^Ri(A',z)(Pi+zi)-^i?i(A',z)dAdA'. 


1 - Ai 1 - Ai 


1-Ai 


(2i7r)2 

We deduce as in the proof of Lemma 5.6 that Vz( 7 zlT 7 z) is trace class, which concludes the 
proof. □ 


iThe proof in [15] is actually stated for operators acting on Lr'(R^). However, the proof applies straightforwardly 
to our bounded domain case Lper(r). 
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5.5 Proof of Proposition 3.1 and Corollary 3.2 

We now proceed with the proof of Proposition 3.1. The assertion (3.3) was proved in Lemma 5.7. 
To get (3.2), we write that 

lT( 7 lt) = ^Iil(-A 7 )+ 1 t (P 7 ) and = , 

so that 

\ ^ |1il(-A7) -ITl (-Al7)I + ■ 

The proof of (3.2) then follows from Lemma 5.6 and (5.24). 

We now prove Corollary 3.2. We compare the total energies 
^peT(7) - ^ (1 t(-A7) - 

T 2 (^1 (^7 Mperi P Pper) D\ Pperi P^l Pper)) j (5.26) 

and notice that 

|l^l(P7 ~ Ppen P7 ~ Pper) ~ 1^1 (P7i ~ Ppen P^l ~ Pper)| = |l^l(P7 ~ PjL-:P'y “f PlL ~ 2 pper)| 

(5.27) 

Using for instance the inequality (recall that |a* I < |a^| < la^l) 


yf,geLl^fT), \Dff,g)\ = 

E 

Ck (/)ck ( 5 ) 

|l,|2 

< ^ V 

~ Ia*|2 

Ck (/)ck (g) 


RG7?,\{0} 

|k| 

1*^31 Rg7?,\{0} 



- I *|2 ipi ll/lli^er(r)ll5llL=„(r)’ (5-28) 

laal hi 

and combining (5.26), (5.27) and (5.28), we obtain 
|^peT(7) - ■S’i'’"(7z,)| < \ |Il(-A7) - 

+ 2|a*|2|r| ~ -f"por(r) ~ 2Mper||L2^7r) • 

Corollary 3.2 is therefore a consequence of Lemma 5.6, (5.24) and the embedding L“j.(r) 
l^^er(r). 

6 Proof for the nonlinear reduced Hartree-Fock case 

In this section, we prove the exponential rate of convergence of the supercell model to the periodic 
model in the nonlinear rHF case (see Theorem 3.3). The proof consists of three steps. 

Step 1: Convergence of the ground-state energy per unit volume 

In the sequel, we denote by Vq := (p^^ — Pper) *rGi and Vl^ '■= {p-y^ o ~ Pper) *rGi (see also (2.6) 
and (2.12)). We recall that 

Hq = ~^A + Vq and 70 = l(i 7 o < ^f) act on L^(IR^), 

HLfi = “^Al + 14,0 and 7 l.o = < £f) act on LI^^{Tl). 
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We denote by 5 > 0 the gap of Hq around the Fermi level ep- 

It was proved in [4] that the sequence converges to Vq in L“j.(r). We will prove 

later that this convergence is actually exponentially fast. As a result, we deduce that for L large 
enough, say L > the operator is gapped around ep, and one may choose the Fermi 

level of the supercell Sp dehned in (2.6) equal to ep. We denote by gp the size of the gap of Hp^ 
around Sp. Without loss of generality we may assume that L^ap jg large enough so that 

VA>AS"P, 5l>|. 

In the last section, we proved that the constants C G K"*" and a > 0 appearing in Proposi¬ 
tion 3.1 are functions of the parameters TZ, |jI4||r,<=a, g and ep of the problem only. In particular, 
it is possible to choose C G K"*" and a > 0 such that, for any choice of potentials V among 
|Vo, (^i,o)r,>igap the inequalities (3.2), (3.3) and (3.4) hold true. 

We first consider V = Vq in Proposition 3.1. We denote by 'yp G Vp the one-body density 
matrix dehned in (5.6) for this choice of potential. Together with Corollary 3.2, we get 

VL e N*, = L-3£:^p-(^r.,o) < ^"^£^'(7^) < fp^eT(7o) + Ce-“^ = + Ce-“^. 

On the other hand, choosing V = Vp^ with L > in Proposition 3.1, and denoting by 7 )^ G Pper 
the one-body density matrix dehned in (5.6) for this choice of potential, we get 

VL > Ls-P, (70) < (7l) < (7l,o) + Ce'^^ = + Ce'^^. 

Combining both inequalities leads to 

VL > LS'^P, |T-3/)(p" - /^p-1 < Ce""-^. (6.1) 

This leads to the claimed rate of convergence for the ground-state energy per unit cell. 

Step 2: Convergence of the ground state density 

In order to compare p^g and p-yj^ g, it is useful to introduce the Hamiltonian := --Ap + Vq 

acting on Lpgj.{Tp). We also introduce jp := 1{H^ < sp). Note that 7 ^ G Vp is the operator 
obtained in (5.6) by taking V = Vq in Proposition 3.1. Therefore, according to this proposition, 
there exist C G and a > 0 such that 

VLGnf \\pyg-pyg\\p^(r)<Ce-^^. ( 6 . 2 ) 

In order to compare py^ with Py^ g, we note that, since 7 l_o is a minimizer of (2.5), then, using (3.4) 
and (6.1), we get that, for any L G N*, 

0 < L-3£:^p-(7i) - L-3£^'’"(7i,o) = {L-^S^^-{^p) - E^irilo)) + ^"( 70 ) - ^-^^^-'(yL.o)) 

< 2 C'e-“-^, 


so that 


VL gN*, 0 < £^"(7^) - ■S’r"(7i.o) < L^2Ce-°‘^ < C"e"“'^ 


for some constants C G K"*" and a' > 0 independent of L. This inequality can be recast into 

VL G N*, 0 < V^p {{Hp^Q-ep) (7L - 7 l.o)) + ^-01(^71. - - Pil,o) < 

Both terms are non-negative, so each one of them is decaying exponentially fast. From the in¬ 
equality (recall that we assumed [a)" I < |a^l < |a5l) 


VfGLl^iT), ||/*rGi||i.^^(r)= ^ 


|ck(/)r 


< 


E 


ick(/)r 


RG7?,\{0} 1^1 l^sl Rg7?.\{0} 


yDiifJ), 
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we obtain that 


VLgN*, WiPjr. - P'fL,o) *r GiWh (T) < T^DiipL-PJ^^„,PL-Pjr..o) < (6-3) 


- |a5'^ 


Consider W G Lpgj.(r). It holds that, for any L > L^ap^ 
J^iPiL - PiL.oW = ^ J2 [((7l)q - (7l,o)q) bh] 

1 


QGAi 


2 mL^ 
1 

2 \ttL^ 


QeA 




(r) 


((Pt-l.o - Pil) *r Gi) 


-W dA 


(6.4) 


E 4 T^L^edr) Q)t^ ((P7.,o - P 7 J *r Gi) B 2 ^(A, Q)-^iy] dA, 

QeAi,”''^ 7 1 1 / 

where is the operator defined in (5.8) for H = Hq, and i?|' is the one for H = Hl q. From 
the expression of the constant Gi in (5.10), we deduce that there exists a constant Gi G M"*" such 
that, for all L > 

VAg-^, VQGAi, ||i 32 P“(A,Q)||g(^,^^(p)) <G) and \\B^iX,Q)\\siLi„(T)) 

As a result, 


J {PjL P7 l,o)^^ 


< 


mGi 

2 tt 


1 


1-Ai 


{(PiL.o - Pil) *r Gi) 


e 2 (i^„(r)) 


1 


1 - A 


-W 


e 2 (i^„(r)) 


We deduce from the Kato-Seiler-Simon inequality [15, Theorem 4.1] and the estimate (6.3) that 
there exists constant G G M+ and a > 0 independent of IF such that. 


< Ge-“^ IIW^IIl^^d ■ 


J {PlL P7i,o)^ 

This being true for all IF G Lpgj.(r), we obtain 

VL>LS’"P, < Ge 


— OlL 


This proves the convergence in Lpg^(r). To get the convergence in (T), we bootstrap the 


procedure. Since (p^^ „ - p^J G Tper(r), then - p^^) *r Gi G T“r(r) with 


VL > ||(p^,,g - p^J *r Gijj^^ (T) < G'e 


per \ 
f ^—aL 


(6.5) 


Consider IF G Lpgj.(r). Performing similar calculations as in (6.4), we get (with obvious notation) 
J^iPiL -PiL.oW 

= ^ E Tr^^„(r)(^r(A,Q)((P7r.o-P7J*rGi)i?|'(A,Q)-^lF-^)dA, 

QeAi •1'^ \ 11 / 

so that 


—2 

/ l^l Cl 

{PiL ~ PiL.o)^ — 11(^71 ,,0 ~ Pil) *r Gi| 


i?°edr) 


^ -IF- ^ 


1 - Ai 1 - Ai 


ei(i" (C) 


and we conclude from (5.25) and (6.5) that there exist constants G G and a > 0 such that 

< Ge""-^. 


VL > LS'^P, 

11 PlL PlL 

Together with (6.2), we finally obtain 


VL > LS'^P, 

11/^70 ~ PlL 
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Step 3: Convergence of the mean-field Hamiltonian 

Finally, since 

Hl ~ Ho — ip^L.o ~ Pio) *r G\, 

the estimate (6.5) implies the convergence of the operator — Hq to 0 in with an 

exponential rate of convergence. 

Remark 6.1. The convergence of the operators implies the convergence of the eigenvalues. More 
specifically, from the min-max principle, we easily deduce that 

sup sup |e„,q[iFL] - en,q[Ho]\ < Ce~°‘^ 

q^r* nGN* 

where denotes the eigenvalues of the operator Hq ranked in increasing order, count¬ 

ing multiplicities. 


7 Numerical simulations 

In this final section, we illustrate our theoretical results with numerical simulations. The simula¬ 
tions were performed using a home-made Python code, run on a 32 core Intel Xeon E5-2667. 


The linear model (Proposition 3.1) 

We consider crystalline silicon in its diamond structure. A qualitatively correct band diagram of 
this system can be obtained from a linear Hamiltonian of the form H = — f A -|- Hp“, where the 
potential Hp™ is the empirical pseudopotential constructed in [6]. The lattice vectors are 

ai = ^(0,1,1)^, a 2 = ^(1,0,1)^ and a 3 = ^(1,1,0)^ 

and the reciprocal lattice vectors are 


at = , 




and 


"fa.i.-iy, 


where the lattice constant is [6] a = 10.245 Bohr (that is about a = 5.43 A). The pseudopotential 
Vp^“ is given by the expression [6] 


where 


l/p'™(x) = ^ t4e‘'^ ’‘ with VkeT?.*, 

kGK* 


Vk = 5'[k] cos 


a{ki 


ko) 


<5[k] 


-0.105 if 

0.02 if 

0.04 if 

0 otherwise. 


|kp = 3(27r/a)^ 
|k|2 = 8(27r/a)2 
jkp = ll(27r/a)^ 


(7.1) 


This system is an insulator when the number of particle (electron-pairs) N per unit cell is iV = 
4, so that the hypotheses of Proposition 3.1 are satisfied. In the sequel, the calculations are 
performed in the planewave basis 


X = 



k G n*, 



^ Acutoff r 5 


(7.2) 


where the cut-off energy is Acutoff = 736 eV. The corresponding size of the basis is |X| = 749. 


In Figure 2, we represent the error on the ground state energy per unit cell and the 
error on the ground state density (in log scale) for different sizes of the regular grid. The value 
of L in (4.7) varies between 4 to 28. The quantities of reference are the ones calculated for the 
regular grid of size 60. We observe in Figure 2 the exponential convergence for both the energy 
per unit cell and the density as predicted in Proposition 3.1. 
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Figure 2: The error on the ground-state energy (in eV) and the L°° error on the ground-state 
density with respect to the size of the regular mesh for the linear model. The logarithm of the 
errors are represented. The linear regression curves are also displayed. 


The rHF model (Theorem 3.3) 

We now consider the rHF model. To our knowledge, no pseudopotential has ever been designed 
for this model. Since constructing pseudopotentials is a formidable task, we limit ourselves to the 
following poor man’s solution, which does not aim at capturing the physics but only at illustrating 
numerically our theoretical convergence results. We decompose the potential self-consistent Vb 
appearing in (2.12) into 

^0 = (P70 ~ Mper) *r Gl = P70 *r Gl — y^per *r Gi, 

and we make the approximation Vq = Fp'™, where Fp“ is the pseudopotential defined in (7.1). 
This leads to the rHF pseudopotential of the form 

:= C - P70 *r Gl. 

In practice, we calculate Fp™ with the potential obtained previously for the grid of size 60. The 
minimization problem (2.5)-(2.6) is solved self-consistently in the basis X defined in (7.2) (we refer 
to [3] for a survey on self-consistent procedures for such problems). We stop the self-consistent 
procedure when the difference between two consecutive densities is less than 10“^. The 

size of the regular mesh varies between 8 to 36. The quantities of reference are the ones calculated 
for the regular mesh of size 60. The error on the energy per unit cell and the error on the 

density are displayed in Figure 3. 

We observe in Figure 3 the exponential convergence announced in Theorem 3.3. 
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Figure 3: The error in energy (in eV) and the L°° error of the density with respect to the size 
of the regular mesh for the rHF model. The logarithm of the errors are represented. The linear 
regression curves are also displayed. 
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